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A CUSS OF UNSTEADY, THREE-DIMENSIONAL FLOW STRUCTURES 
IN TURBULENT BOUNDARY UYERS 

By 

Robert L. Ash^ 

SUM4ARY 

A restricted class of mathematically admissible, imsteady, three- 
dimensional flows has been identified which may constitute part of the 
structure observed in turbulent boundary layers. This report presents 
the development of the model and discusses some general results. The 
resulting solution has characteristics vdiich suggest how upwelling low- 
speed flow can trigger a downward jetting of irrotational high-speed 
fluid into the boundary layer. 
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unit vectors 
pressure 

local velocity vector 

tree-stream velocity 

v«)locity components 

complex velocity function 

reciprocal streamwise disturbance length 

ratio of disturbance to free-stream velocity 

kinematic viscosity 

complex variable 
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INTRODUCTION 


Experimental investigations continue to show that turbulent boundary- layer 
flows are composed of coherent, unsteady flow structures distributed randomly 
in space and tiM (see ref. 1, for example). The mechanisms responsible for 
the evolution and destruction of those structxxres are not understood. Neither 
is the process responsible for the production and distribution of Reynolds 
stresses. In order to modify turbulent flows in a favorable or controlled 
manner, the physics of the flow must be better understood. 

Presently, experimental flow visualization studies have shown that large, 
unsteady vorticular structures are a major part of the outer turbulent 
boundary layer (refs. 1 and 2). The actual fora and size ot a single, identi- 
fiable flow structure is unknown because there is no method available to 
separate one structure from the remnants of all the others in which it is 
immersed. Furthermore, Hama (ref. 3) has shown how visual observations of 
unsteady streaklines can lead to erroneous conclusions and to inproper isolation 
and interpretation of structural features. The present investigation has 
attempted to identify mathematically allowable solutions to the unsteady-three- 
dimens:' onal Navier- Stokes equations which may constitute part of a turbulent 
boundary layer. The ultimate objective of this work is to enable experimental 
identification of turbulent flow structures and develop techniques for 
modifying turbulent flows favorably. 

Flow structures —identified loosely as vortices— are apparently embedded 
in turbulent flows. A variety of structures may be present, and none of them 
are presumed so large that they are required to satisfy inner and outer 
boundary conditions simultaneously. Rather, structures near the wall are 
presumed to satisfy wall boundary conditions while outer structures must 
satisfy the irrotational outer conditions. I'he approach employed here was 
to attempt to find solutions to the constant property Navier-Stokes equations 
which could be subjected to either the inner or outer boundary conditions 
(but not both). Thus far, models appropriate to outer flow conditions appear 
to show more promise than the near wall structures. The results of the outer 
stiructure studies follow. 


GOVERNING EQUATIONS 


The local velocity vector u, given by 

u ■ u ej ♦ V e 2 ♦ w e 3 (1) 

where ej, 62 , and 03 are unit vectors in the x-, y- and z- directions of 
a Eulerian coordinate, respectively, tdiich must satisfy the Navier-Stokes 
equations: 

Continuity; 

V • u ■ 0 ( 2 ) 

and Conservation of Momentum: 

~ ♦ (u • ^) u ■ - - * vV^u (3) 

db 0 


where 


^ = 


ei 




(4) 


and P and p are pressure and density, respectively. 

Referring to figure 1, a coordinate system is chosen which has its 
origin at the nominal turbulent -nonturbulent interface. The boundary is 
assumed to be initially parallel, and if no interfacial motion were present, 
the velocity distribution for positive values of y would simply be 


u • U ei 

m 4 


(5) 


In reality, the turbulent -nonturbulent interface is characterized by 
three-dimensional peaks and valleys which occur rather abruptly, then are 
obliterated slowly as they convect downstream (see ref. 1). The mathematical 
model which follows has attempted to couple the low-speed burst disturbances 
from the wall region with outer eddy structures. It is assumed that the low- 
speed fluid parcel ejected from the wall beneath a flat plate, turbulent 
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boundary layer is not able to penetrate the turbulent-lnviscid interface. 
Rather, this slow moving fluid is sxirrounded by fine structure turbulence; 
as it arrives at the upper reaches of the boimdary layer, the fine structure 
shields the streaky fluid from the free stream. The large velocity gradients 
which anist be present act to sustain the fine structure. In addition, the 
upward moving filament causes a small bulge to occur in the turbulent -non- 
turbulent interface region. 

Consideration of pressure changes resulting from a bulge suggests from 
Bernoulli's equation that, if the velocity above the disttirbance increases, 
the pressure decreases, drawing the bulge further upward. This type of 
disturbance would tend to pull the turbulent fluid rapidly upward into the 
irrotational regime as it is swept downstream. However, there appears to 
be no stopping mechanism to control bulge height. An alternate model would 
be to assume the inviscid fluid maintains a constant volume flow by pro- 
ducing an interfacial depression to compensate for the bulge. Under those 
conditions, the bulge need not produce a local velocity increase and the 
destabilizing pressure disturbance need not occur. Those disturbances have 
been studied here. 

Consider the class of velocity distributions given by 
u - ei ♦ v(x,y,2,t) %2 * w(x,y,z,t) 03 (6) 

where is assumed constant. Then, the continuity equations reduce to 


9v ^ 3w 
3y 3z 


0 


(7) 


vdiich can be satisfied when v and w are functions of x, t and complex 
variable ? , where 


C ■ y ♦ i z 


( 8 ) 


Then, eqxiation (7) becomes 


3v . 3w 
3 ^ " 3 ? 


( 9 ) 
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which can be integrated to yield 


w(x,5,t) ■ i v(x,C,t) f g(x,t) 


( 10 ) 


Here, we choose g(x,t) ■ 0, since no apparent benefits result from: inclusion 
of a function of x and t alone. 

Equation (10) can be utilized in equations (3) and (6) to develop ex- 


pressions 

for 

the conservation 

of I 

X or 


direction: 
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■ n 




3X 

■ V 
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62 < 

direction: 
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A !1 
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32v 

3t 

♦ U 

■ 3x " ’ p 

V 

35 

3x2 

and 3 

: or 

03 direction: 



( 11 ) 


( 12 ) 


ifil 3v . .ilP ,, ifv-l 
^ Lat « 3x p 35 3^2 J 


(13) 


Obviously, equations (12) and (13) are the same. Furthermore, the re- 
quirement that P- ■ 0 implies that p • p(5,t), which requires that pressure 

0A 

fluctuations consistent with the above equations be "global" in the x- 
direction. Those pressxire fluctuations are not particularly interesting, and 

it is therefore assumed that • 0. Hence, the equations of motion are 
satisfied by 


3t « 3x 


3fv 

3x2 


(14) 


where the nonlinear terms in the momentum equations are cancelled through the 
relationship: 


iv(x,5,t) 


(15) 
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The linear partial differential equation (14) can be solved using clas- 
sical techniques, and if V is defined by 

V • V - iw (16) 


the solution to equation (14) subject to initial conditions: 

V(x,it,o) ■ v(x,i:,o) - iw(x,i 2 ,o) (17) 


can be satisfied by 

J*" *xp v(7.ic.o) «xp (i(«c»en)ldTd^B<Ui (18) 

where the overbars on the spatial variables have been made dimensionless by 

U.( ) 

( ) . (19) 


and dimensionless time is given by 


t 



( 20 ) 


Presently, the initial velocity disturbances which are considered global 
through equation (18) are not as interesting as spatially periodic. That is, 
if the disturbance is finite with respect to the spanwise coordinate, the 
velocity distribution is given by 


v » f (x,t) exp (i:->:] 

and 

w ■ if(x,t) exp (i:-y)] 


(21) 


( 22 ) 


where is an unspecified spanwise wavelength, and f(x,t) is a solution to 



•t 'X ix- 
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with 



( 24 ) 


The solution to F can be gotten using Fourier transforms as 


F(x,t) 




e^®^ F(C.o) dedB 


(25) 


At this point, solutions to F(x,T) are being studied for various types 
of initial disturbances in the x-direction. The results of that investi- 
gation will require considerably more computational time and development of 
additional computer graphics capabilities in order to verify the presence of 
these structures in a turbulent boundary layer. 

One disturbance which has produced interesting results is 

F(x,o) - r e'^^‘ (26) 


There, the integration of equation (25) is possible and 


i- . '.L i-. i/ril'-il cS-^-i^s) (27i 

(4y2?*1)1/: 

Under those initial conditions (a small area of fluid is given an initial 
upward velocity distribution at the turbulent-nonturbulent interface), the 
model predicts that the bulge is accompanied by a laterally shifted downward 
jetting of high-speed irrotational fluid into the boundary layer. This 
process may approximate the inrush phase near the wall, which triggers new 
bursts. 


CONCLUSIONS 

A mathematical model has been developed which predicts deformation of 
the interface between a turbulent boundary layer and the outer, irrotational 
flow. Although solutions can be developed which describe the full x-y plane 
initial velocity conditions, those solutions do not appear to yield results 


f 

,4T 

I 

► 

! 

which describe physics within the turbulent boundary layer. On the other 
hand, a variety of restricted solutions can be developed whicn represent 
surface deformations resulting from initial velocity disturbances in finite 
regions of the x-y plane. Those calculations have shown that bulges at 
the interface can produce jetting of irrotational fluid toward the wall which 
may describe some turbulent features. Investigation of the finite area 
disturbances should be pursued to permit either isolation or elimination of 
the proposed structures a.s important features in a turbulent boundary layer. 
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